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Abstract 

The physical foaming agent is dissolved in the polymer at a pressure and temperature beyond its critical condition and the foam 
structure is formed by a sudden pressure drop or a temperature increase in the mixture in three steps; nucleation, growth and 
coalescence. In the nucleation step, due to a supersaturation, thermodynamic unstability is formed in the mixture. This induces 
a tendency in the solvent molecules for a phase transition from a supersaturated state to a gas state in the direction of 
unstability reduction. By overcoming the energy barrier, the free energy of the system is reduced and gas stable nucleuses are 
formed in the free volumes of the polymer chains. Growth and coalescence of the nucleuses take place with the diffusion of the 
gas molecules inside them. With increase in the nucleation efficiency, the growth and coalescence rate are reduced. Therefore, 
prediction of nucleation rate is an important factor for controlling the thermoplastic foam structures. Investigation of the 
nucleation step has been doable with the aid of nucleation theories. The classical nucleation theory is a prominent method for 
the investigation of nucleation phenomena in thermoplastic foams, but due to the existing of a divergence between its 
theoretical values with the experimental data, some modifications have been introduced into this theory. Other competing 
theories such as the Density Functional Theory and the Self Consistent Field Theory have also been developed. The main aim of 
this review paper is investigation of the existing theories and their modifications for thermoplastic foams. 
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Introduction 

Foams are materials in which gas cells are surrounded by a dense matrix. The first generation of thermodynamic 
foams had lower ductility and strength compared to pieces of solid polymers. This deficiency was resolved by 
reducing cell sizes in foam structures. In fact, cells with smaller dimensions in relation to a critical size can inhibit 
crack growth in the bulk. Thus, by developing such cells inside the material it is possible to reach pieces with 
decent mechanical properties and less material consumption [1, 2]. Based on this fact, investigation of the effective 
parameters for controlling cell size in thermoplastic foams has been carried out extensively and resulted in the 
production of microcellular foams. There has also been a considerable tendency for the development and 
production of nanocellular foams [2-18]. 

Both physical and chemical foaming agents can be used in the two common methods for the production of 
polymeric foams. When chemical foaming agents are incorporated, the foam's structure is usually formed in the 
polymer bulk due to a chemical reaction. In the first method such as polyuretan foams. Chemical foaming agents 
are the byproduct of the reaction between the polymer constituents in the form of gas which are trapped between 
the polymeric chains. In another method, such as aerogels and thermoset foams, the foam structure is formed by a 
chemical reaction resulting in a crosslinked three dimensional network [19-22]. On the other hand, when physical 
foaming agents are incorporated, it is dissoleved in a polymer at a pressure and temperature beyond its critical 
condition and the foam structure is formed by a sudden pressure drop or a temperature increase in the mixture 
[23-25]. 

Physical foaming agents are used in continuous and batch processes in extruders and pressured vessels 
respectively, for the production of thermoplastic foams [2]. 

The foaming process is carried out during the following steps: nucleation, growth, coalescence and structure 
stabilization. After nucleation of gas molecules and development of stable nucleus, gas molecules diffuse inside the 


11 



www. seipub . org/papt 


Plastic and Polymer Technology (PAPT) Volume 4, 2016 


nucleus and make them grow into cells. In case of a continuous growth, cells become connected and result in 
bigger cells [26, 27]. 

The main goal in this research is to review and discuss the current theories for the prediction of nucleation rate in 
thermoplastic foams produced by physical foaming agent by a sudden pressure drop. 

The final foam structure is dependent on factors such as size, density, size distribution, dispersion and coalescence 
of cells. With increase in the efficiency of the nucleation stage, an additional number of molecules become involved 
in the nucleus production, thus on one hand there is less chance of cell growth and coalescence and on the other 
hand, there is an increase in cell density and a decrease in cell size. Therefore, this stage is a controlling factor in the 
foam structure. Study of the nucleation theories and recognitions of their strengths and weaknesses can 
significantly assist the control and prediction of the thermoplastic foams structure [28]. 

Mechanisms of nucleation phenomenon in thermoplastic foams are investigated using nucleation theories that are 
used to describe nucleation in first order phase transitions such as crystallization and boiling. 

With the occurrence of super saturation and chemical potential difference in these processes, molecules are locally 
arranged beside one and other and fluctuations of energy and density are raised in the system. In case of exceeding 
the energy barrier, the free energy of the system is reduced and stable sites are formed in the first phase from the 
new phase which is followed by nucleation and phase transition. In thermodynamic science this phenomena is 
mainly modeled by classical and non classical methods. The aim of these theories is to investigate the effects of 
thermodynamically and kinematical variables and physical and chemical properties of constitutes on the 
nucleation. The nucleation rate is predicted through the calculation of critical nucleation energy and the critical 
sizes of the newly formed regions in the new phase [29]. 

In the classical nucleation theory the new phase is assumed as a nucleus with a critical size which is separated 
distinctively from the former phase by a dividing surface. 

In other words, the system is supposed to resemble an ideal model consisting of two homogeneous phases; old and 
new, and an imaginary dividing surface at the border between the phases. The thermodynamic properties of this 
system are a result of the thermodynamic properties of the two phases and the dividing surface. Additionally, the 
free energy of the system includes a volumetric part (two homogeneous phases) and a surface part (the dividing 
surface) [22, 23]. 

This theory consists of both homogeneous and heterogeneous nucleation. Homogeneous and heterogeneous 
nucleation occurs in the single phase primary mixture and in the presence of a second phase consisting of additive 
particles respectively. The second phase plays the role of the nucleation sites which cause a reduction in the 
nucleation energy and an enhancement in the nucleation rate [2, 29- 35]. 

The nucleation theory is applicable commonly for the calculation of nucleation rate in thermoplastic foams in the 
presence of a supercritical solvent. Although, the theoretical and experimental results do not show a proper 
consistency and in most cases only a correct pattern of changes may be predictable. Scientists have declared two 
main reasons for this inconsistency; error in the calculation of variables and simplifying assumptions for the 
establishment of the relations between variables. The variables include the interfacial tension between polymer and 
solvent, difference between the inside and outside pressure of the nucleolus. Other sources of error in this theory 
which should be taken into account are: absence of the effect of time, pressure drop rate and super saturation 
content. A lot of modifications have been suggested by various scientists to reduce the inconsistency. These 
attempts have been in a direction to care for the mentioned sources of errors which in some cases have resulted in 
superior outcomes [3, 36-41]. 

The presence of a number of small molecules in the nucleus and the assumption of an existing dividing surface on 
the border between the two phases raises uncertainty which cannot be resolved by this theory. To begin with, what 
exactly is a cluster? Are there any criterions physically for positioning the dividing surface? How can one 
determine the actual shape, structure and density of small clusters? How accurate is the use of specific surface 
energy as an independent factor from cluster size [29]? 

Therefore in order to eradicate the mentioned misunderstandings in the field of thermodynamics, some 
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competitive theories with different approaches for the interpretation of first order phase transition have also been 
used from which we can point out to the density functional theory (DFT) and self consistent field theory (SCFT). 
The density functional theory is the basis of non classical methods for modeling of nucleation in the first order 
phase transition. In this theory, nucleation is considered as a local fluctuation of molecules' density. The free 
energy of the system is calculated on the basis of the fluctuations of density or weight fraction of the mixture 
components through the elimination of the dividing surface between the old and the new phases [38, 42-46]. It 
should be mentioned that this theory is used more often in the case of small molecules, but recently there have 
been some reports of its application in polymeric solutions for the calculation of mixture free energy and interfacial 
tension of solvent and polymer [45, 46]. 

In the SCFT theory, change in free energy is calculated versus the weight fraction of the mixture components and is 
used for the modeling of nucleation in polymeric foams. This theory is a mean field theory and is established for 
the study of non-uniform polymeric systems based on the equilibrium statistical mechanics in which there is no 
mention of time. In this theory, potential energy models such as Lenard-Jones (L-J) and mobility models such as 
Gaussian Chain model have been applied for the effects of molecular interaction between solvent and polymer and 
chain movements, respectively which will finally result in the derivation of free energy equation for a bubble in a 
compressible polymeric fluid versus the weight fraction. Due to the fact that these theories are based on molecular 
models, their findings have only been used for studying the effect of parameters, deeper understanding of 
nucleation and efficiency enhancement and thus are not comparable with experimental results [3, 47-53]. 

Nucleation in thermoplastic foams has also been studied from the perspective of phase separation kinetic in 
compressible polymeric solutions. Phase separation in polymeric solutions occurs with deviations of temperature, 
pressure and weight fraction of mixture components from the equilibrium state. The homogeneous mixture enters 
the unstable and metastable regions during phase separation and this phase separation takes place under 
nucleation, growth and spinodal decomposition mechanisms. This route depends on the amount and how of 
changes of temperature, pressure and weigth fraction of mixture components. In a case where the system is pushed 
towards the totally unstable thermodynamic region the phase separation is the spinodal decomposition kind. If the 
system is pushed towards the metastable thermodynamic region which is between the binodal and spinodal lines, 
the phase separation mechanism is in the form of nucleation and growth. The spinodal decomposition results in 
the formation of a bi-continuous or a percolation microstructure. Meanwhile, the formation of separated regions in 
the matrix is due to the mechanisms of nucleation and growth [54-59]. 

Some models have been developed for further investigation of the behavior and kinetic of phase separation of 
compressible polymeric solutions. Regarding the interactions between gas and polymer molecules (e.g. formation 
of hydrogen bonds) in polymeric solutions, the mentioned models that include such effects is considered recently. 
The difficulties of the calculation of experimental parameters which are used in the models for the prediction of 
phase behavior of compressible polymeric solutions; application of models with less experimental parameters is 
more convenient [60-66]. 

There are also some reports on the use of light scattering and neutron scattering techniques for the study of the 
kinetics of phase separation in compressible polymeric solutions [67-69]. 

Due to the broadness of nucleation theories in the field of thermoplastic foams, the context of nucleation in 
thermoplastic foams from the point of view of phase separation kinetics in compressible polymeric solutions will 
not be reviewed. Thus, here an attempt has been made to initially study the nucleation phenomena and classical 
nucleation theory in thermodynamics and in the next step their applications, limitations and modifications for 
thermoplastic foams has been covered. Finally the density functional theory and Self Consistent Field Theory with 
their applications in polymeric solutions and theormoplastic foams have been reviewed briefly. 

Nucleation in First Order Phase Transitions 

Nucleation is the primary stage of first order phase transitions such as crystallization and boiling in which a 
substance can be transformed in to a stable state from a semistable. The transition from a metastable to an actual 
stable state occurs in the fluid to have a lower-enrgy state [22]. Generally, the thermodynamic deriving force for the 
first order phase transition, chemical potential difference or supersaturation, is 
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Al* (G 0 id G new )/M P old P new (1) 

Supersaturation (Ap) is the gain in the free energy per molecule or (atom) in refeternce to the passage of the phase 
from the minimum with higher Gibbs free energy (G 0 id) to the minimum with lower Gibbs free energy (G new ). g 0 i d 
and g new are, the chemical potentials of the old and new phases at the crospending minima respectively [29,30]. 

For nucleation and phase transition of material in a specific period of time, a state of supersaturation as well as an 
energy level higher than the energy barrier is required [29]. 

As an example, the transition of the Van der Waals (VDW) fluid with below equation of state from liquid phase 
(old phase) to gas phase (new phase) for better understanding of first order phase transition is expressed. 

(P + M 2 a'/V 2 )(V - Mb') = MkT (2) 

Where P and V are fluid pressure and volume, respectively, M is the number of molecules (or atomes) in this 
volume, T is the absolute temperature, K is the Boltzman constant, and a and b are material constants acconting, 
respectively, for the molecular interactions and molecular volume. These two constants can be written in terms of 
the critical pressure P cr volume V cr and tempreture T cr of the VDW fluid using the following relations: 

P cr = a'/27b' 2 V cr = 3Mb' T cr = 8a’/27kb' (3) 

The VDW fluid equation of state and free Gibbs energy can be expressed by equations 4 and 5, respectively, using 
the appropriate relations for reduced pressure P = P/P cr , volume V = V/V cr and temperature T = T/T cr . 

(P + 3/V' 2 )(3V'- l) = 8T (4) 

G(V ,T’) = G ref (T ) + MkT cr [3 V T /8 — 9/8V — T ln(3V — l) (5) 

Nucleation in phase transition from a liquid state to a gas state in constant temperature with volume changes has 
been investigated using the phase diagram of the VDW fluid and changes in energy shown in figures 1 and 2. 



FIG 1 . DEPENDENCE OF REDUCED PRESSURE OF VDW FLUID FIG 2. DEPENDENCE OF GIBBS FREE ENERGY OF VDW FLUID 

ON REDUCED VOLUME [29] ON REDUCED VOLUME [29 

In case the initial material phase is in the state of liquid, with a certain temperature and pressure and a volume 
equivalent to point a in figure 1, its volume has to increase to point h in order to undergo a phase transition to the 
gas phase in the same temperature and pressure. The most probable phase transition rout is through a uniform 
change in material density and an increase in its volume. According to figure 2, the free energy increase initially, 
reaches a maximum and reduces in the route with an increased volume. So the free energy difference of points a 
and maxima is the energy barrier of uniformly changing density route. Considering the amount of required 
energy for exceeding the energy barrier of the mentioned route, it is understood that phase transition is 
thermodynamically unfavorable [29, 30]. 

Another conceivable phase transition path with much less energy consumption in comparison to the first path is 
known as the non-uniform change of the density of the old phase into the density of new one through local 
changes and fluctuation of density in a contracted spatial region by a small number of molecules. 
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If M and n*are, respectively, the number of total molecules and the number of molecules in the spatial region of 
locally changed density, the energy barriar for the phase transition corresponding to the nonuniformly changing 
density (n* Ap) is much less than the energy barriar for the path corresponding to uniformly changing density (MAg) 
and thus the possibility of first order phase transitions to occure by local density fluctuations is generally greater 
than by uniform change of the density of the old phase. 

Many observations such as rains, snowfalls and experimental investigations suggest that the path followed in 
reality by a first-order phase transition is infect the appearance of nanoscopical formations with density close to 
that of the new phase which are randomly formed in the old phase and the termination of phase transition by the 
overgrowth of the new phase to larger sizes. So, what is meant by nucleation is actually the process of random 
generation of those nanoscopically small formations of the new phase which are capable of irreversible overgrowth 
to macroscopic sizes. 

As stated previously, mechanisms of nucleation phenomena in thermoplastic foams have been studied from a view 
point of first order phase transition of solvent molecules (gas) from a supersaturated state to a gas state in a 
polymeric matrix. The existing theories have been used effectively in this field. In continue, first the classical 
nucleation theory and then its application in thermoplastic foams will be looked. 

Classical Nucleation Theory (CNT) 

In the Classical nucleation theory, nucleation presumably occurs in the new phase in the form of nucleus with a 
critical size and it is separated from the old phase by a dividing surface and thus the molecules of the old and new 
phases are separated. 

In other words, after a phase transition the system is supposed to be made up of two uniform and distinguishable 
old and new phases and a dividing surface as demonstrated in the model in figure 3. 



FIG 3. SCHEMATIC VIEW OF OLD AND NEW PHASES AND DIVIDING SURFACE 

In this theory, two types of homogeneous and heterogeneous nucleation exists which are distinguishable. 
Homogeneous nucleation takes place in the initial single phase mixture and heterogeneous nucleation also occurs 
in the presence of a second phase such as additives. The second phase plays the role of nucleation sites and reduces 
the nucleation energy while increasing the rate of nucleation [30]. 

In this theory, the mechanism of changes in the system free energy due to the formation of a nucleus from the new 
phase in the old phase is investigated in equilibrium thermodynamics. Based on this theory, in the case of a 
uniform liquid serving as the initial phase, due to phase transition it will be transformed into a nucleus with a 
distinct boundary in the gas phase which surrounded by the initial liquid phase. 

In general, the total energy of a homogeneous part of the system is a sum of the energies of the surface (no volume 
V = 0), and bulk (no surface area A = 0) contributions. 

U = T 0 S-PV + aA + 'Zn i N i (6) 

Which U , T 0 , P , S , V , A , a, Hi , N t are internal energy, temperature, pressure, entropy, volume, surface area. 
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interfacial tension, chemical potential and the number of molecules of system components respectively. 

So the initial energy of the homogeneous liquid is 

U 0 = T 0 S 0 -P 0 V 0 +? l nl l N l °i (7) 

Which U 0 , T 0 , P 0 , V 0 , [i l t , N i t are the thermodynamic characteristics of homogenous liquid phase. 

The final energy in the system consisting of the vapor cluster surrounded by liquid is 

U = T 0 S - PgV g - P L V L + a A + Z p ig N ig + Z n iyl N tl + Z p is N is (8) 

Which P g ,V g , Hig , N ig , er are pressure, volume, chemical potential , number of molecules of cluster and interfacial 
tension between liquid phase and cluster. V L , P L are volume and pressure of liquid phase and n is , N i s are chemical 
potential and number of molecules of surface respectively. 

When a droplet is formed in a vapour, three choices exist for the free energy which depends on the conditions of 
the real nucleation process (in nature or laboratory) or the hypothetical mind experiment. In turn, the suitable free 
energy is related to the quantities which are held constant when the droplet forms. 

If the pressure of the liquid is constant P l — P 0 and the total molecular numbers are unchanged N i tot = N i L + N i g + 
N is = N-\ then the Gibbs free energies in the initial and final states are 

G 0 = U () ~ TgSg + PgVg = £ fa ° (9) 

G ~ Ug- TgS + PgV tot = U ~ TgS + Pg (Vg + V L ) = (Pg ~ Pg)Vg + A + £ Pi.g + Z Pl,l N iJL + Z Pi.S (10) 

The Gibbs free energy change in the droplet formation for equilibrium thermodynamic is 

AG = ( Pg - Pg)Vg + aA + Z m,g Ngg + z H UI N u + z Pi, S N llS - Z Pi, i Ngi (H) 

Regarding the difficulty of measuring the thermodynamic properties of surface and nucleus, a few assumptions are 
set to ease the calculation of nucleation free energy and to estimate the free energy based on the available 
thermodynamical properties which will be reviewed in the following section. 

By assuming that the composition of the liquid phase is constant, which according to the Maxwell equation and 
state of constant pressure denotes the constancy of chemical potential // ; ; t and knowing that the total 
molecular numbers are unchanged, we reach: 

AG = (1 Pg - Pg)V g + ad + z {Pi,g - Ph) N i: g + z (Pi,s - Pi,i) Ni, s (12) 

When the liquid volume is large and the cluster is small, its formation does not significantly change the state of the 
liquid. In this case, equilibrium can be accomplished by setting the derivatives of AG to zero. The equilibrium is 
reached by maintaining the pressure p and chemical potential /v ; at constant levels and considering the 
deviations in accordance to: 

1. V, keeping N i g and N i s constant. 

2. N ig , keeping Nj g ,j i and IV; jS vi constant. 

3. N is , keeping Nj S , j i and N ig v / constant, 
and setting these derivatives equal to zero. 


p* p — ( dA / f 2a A 

F a F 0-{ / dVg )° ~{ r .) 

(13) 

fti.l ~ fti.s ~ fti,g 

(14) 


Equilibrium is on denoted by (*). 

By replacing equations 13 and 14 into equation 12, the change in free energy due to the formation of a nucleus from 
the gas phase in the initial liquid phase at equilibrium is achieved. If the nucleus is assumed to be spherical then 
we have: 
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?/t* ?rr* 4 T 2 4 2 

— - 14, + cr = — - nr + 47T r <7 = - 7rcr r 


AG* = (P 0 -P;)iy + ff M*= r . r . 3- .... - 3 

If, the homogeneous relations of equilibrium radius and equilibrium free energy are[29,30]: 


r = 


A P* 


AG* = 


167 T(J 3 
A P * 2 


(15) 


(16) 

(17) 


Classical Nucleation Theory in Thermoplastic Toams 

As stated previously, in the case of thermoplastic foams, the physical foaming agent (usually a gas) is dissolved in 
polymer at a pressure and temperature higher than the critical condition. Foam is developed through a pressure 
drop of the solution. Incorporation of the classical nucleation theory for the prediction of the nucleation rate in 
microcellular foams was initially carried out by Martini, Youn and Colton in MIT University in the 80s and since 
then it has been used continuously in the field of thermoplastic foams [2, 3, 26, 27, 70]. 

(1) Homogenous Nucleation 

In order to apply this theory for thermoplastic foams, after a pressure drop the resultant mixture is assumed to 
consist of solvent molecules in gas state surrounded by a polymeric layer as the old phase. A Change in free 
energy of the equilibrium state for the formation of spherical gas nucleus in polymer is expressed by equations 
18 and 19 which is similar to relation 14. 

AG = -V AP + Act (18) 

AG = -nr 3 A P + 4nr 2 a (19) 

3 

AP is the difference in pressure between inside and outside of a cluster. This term can usually be replaced by 
pressure difference between saturation and atmospheric pressure in very fast depressurization, r is the cluster 
radius and y is the microscopic interfacial tension of polymer and gas cluster which for a flat surface is usually 
assumed to be equal to the macroscopic interfacial tension between polymer and gas solvent. 

Changes of droplet formation free energy with cluster radius is demonstrated in figure 4. It is realized that the 
interfacial term increases with r 2 and the potential difference decreases with r 3 , thus, the creation of small 
bubbles lead to a free energy increase to the extent where the radius of the bubble nucleated is larger than a 
critical size relative to a maximum free energy. The critical size of cluster and maximum free energy are 
represented by crtitcal radius and critical nucleation energy respectively [1]. 



FIG 4. CHANGES OF NUCLEUS FORMATION FREE ENERGY WITH ITS RADIUS 

Critical radius and energy in equilibrium state are expressed with equations 15 and 16 which are equivalent to 
the extremum point of equation 18. 

Homogeneous nucleation rate on the basis of Boltzmann relation is expressed 

/ = f 0 C 0 exp(-AG* hom /KT) (20) 

Which f 0 and c 0 are frequency factor as a probability of gas molecules interaction to nucleuses and solvent 
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concentration respectively [2, 26, 27, 70]. In nucleation literature, nucleus formation kinetic is investigated on 
the basis of probability of attachment and detachment of molecules to nucleuses as an effect of interaction [30]. 

(2) Heterogeneous Nucleation 

Addition of a second insoluble phase (solid particles) to the solution leads to a heterogeneous nucleation. The 
interfacial tension is reduced and the nucleus develops in the interface with the addition of the second phase 
known as the nucleation agent or nucleation site. As illustrated in figure 5, the nucleus set a wetting angle 0 
with the surface of particle. The relation for stress equilibrium is presented as: 

Vis =Y gs + cos0 Y g i (21) 

y ,y , y are interfacial tension of solid-liquid, gas-solid and gas-liquid phases. Additive particles, solvent and 

AS gS gl 

polymer are formed solid, gas and liquid phases respectively. 

From relation 21 

c„sO= (Y '>“ Y -V y (22) 

/ tgl 


Heterogeneous critical nucleation energy in equilibrium state is 


AG het - -VAp + A gl y g[ + A gl y gs - A LS y LS (23) 

AGLt = AGh om S(0) (24) 

5(0) = Q) (2 + cos0)(l — cos0) 2 (25) 



FIG 5. INTERFACIAL TENSIONS OF HETROGENUS NUCLEATION 


FIG 6. CRITICAL HOMOGENUES AND HETROGENUS 
NUCLEATION ENERGIES 


Because S(0) is always smaller than one, according of equation 24, hetrogenous critical nucleation is smaller 
than homogenous state and its nucleation rate is increasd. Herogenous nucleation rate is 


Jhet = fiCjexp 


(- ac; “/ kt ) 


(26) 


Which h and c x are frequency factor and solvent concentration respectively [2,70 -87]. 


Weaknesses of the Classical Nucleation Theory 

Application of the classical nucleation theory for the assessment of thermoplastic foams nucleation rate in the 
presence of a supercritical solvent in the batch process is common. Nevertheless, the theoretical results do not have 
appropriate consistency with the experimental results and in many cases this theory is only able to give a decent 
prediction about the trends of changes. The reasons behind this inconsistency will be discussed in this section [3-9, 
26-29,38-41, 50]. 
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As explained in section 2-1, the relations applied for the calculation of nucleation rate is based on equilibrium 
thermodynamic of the system, however, nucleation in fact is a nonequlibrium and irreversible process and its path 
is very effective and determining in the final efficiency. Due to the complications of performing the nonequilibrium 
thermodynamic method, scares number of studies has been conducted in this field so far mainly in relation to 
small molecules [88]. 

Scientists identify other factors for this inconsistency which originate from the assumptions for simplification of the 
relations of classical nucleation theory. As states in former sections, the assumptions are as follows: 

1. The thermodynamic properties of nucleus in nanoscale are equivalent to thermodynamic properties of the 
bulk in macroscopic scale (equation 13). 

2. The interfacial tension difference of nucleus in a nanoscale with a curved surface is equivalent to the 
macroscopic interfacial tension difference between two bulks with a smooth surface which is also known as 
the capillary effect. 

3. With assuming an instantaneous pressure drop, the pressure difference between the inside and outside of 
the nucleus is considered to be equivalent to the pressure difference between the initial and the end of the 
process [3, 38, 53, 89-97]. 

The effect of the stated assumptions is even more apparent with decrease in cell size which extends the 
inconsistency between the theoretical and experimental results [3, 38]. Ignoring the factors related to pressure drop 
rate or time, supersaturation level and the presence of polymeric chains in these relations results in the prevalence 
of errors in calculations [40,55,87]. With regard to the mentioned statements, some scientists have made some 
modifications on this theory which will be pointed out in the upcoming section. 

Modifications for the Classical Nucleation Theory 

The source on inconsistency between the experimental and the theoretical results is recognized as the incorrect 
calculation of energy and the nucleation rate for the modeling of nucleation phenomenon with the use of 
nucleation energy relation (equation 17). The incorrect calculation of nucleation energy can be due to two reasons; 
low precision in the models and experiments for the calculation of the nucleation energy experimental parameters 
and not including all the effective factors in it. Thus, the modifications performed on this theory for the 
thermoplastic foams are also categorized into two methods; the exact calculation of variables and inclusion of other 
variables affecting the amount of nucleation energy to the mentioned relation. The modifications on this theory are 
developed extensively in small molecules and there has been an attempt to eliminate the simplifying assumptions 
for the derivation of nucleation energy relation [90-97]. 

For the exact calculation of CNT variables, some equation states such as Non Random Hydrogen Bonding (NRHB) 
and Sanchez-Lacombe (SL) are used for calculation of the nucleation energy experimental parameters [98-99]. 

A tremendous amount of work has been carried out by Colton and Suh, Han and Han, and Goel and Beckman. 
Shafi and Flumerfelt, Ohshima and his co-workers, and Sirupurapu et al work on predicting bubble nucleation 
in thermoplastic polymers using classical nucleation theory [87]. 

Colton and Suh applied the classical nucleation theory to assess the rate of nucleation in microcellular polystyrene 
foams, carbon dioxide gas and zinc stearate particles and became aware of the inconsistency between the 
theoretical and experimental results due to error in the evaluation of interfacial tension. Then, they introduced the 
effects of polymeric chains in the nucleation energy through specific assumptions. 

Colton and Suh also investigated the effect of the presence of soluble and insoluble additives in a solution on the 
nucleation energy. They assumed that, a dissolved additive increases the free volume and distance between 
polymer chains. If the solubility limit of the additive is exceeded, the additive will precipitate and form particles as 
a second phase and nucleation will be heterogeneous. The ability of the particles to act as nucleation sites can be 
calculated by determining the degree of bonding between the particles and the polymer. If there is a high degree of 
bonding, a large amount of energy will be required to force the interface apart. If this energy is high enough, 
nucleation may be precluded. If a poor bond exists, the interface will have a large amount of energy associated 
with it and will enhance nucleation [27,70]. However, they did not take into account the effects of interaction 
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between particles and polymer on the distribution and dispersion quality of particles and its effect on the final 
foam structure. 


Nucleation is affected from local change of free volume with pressure and additive particles in their model. The 
Lennard-Jones model which approximates the potential energy of two non-interacting molecule for modeling of 
interaction energies has been used. The interaction between the electric fields of the negative electrons and the 
positive nuclei give rise to forces of interaction between any pair of molecules, which depend on the nature of the 
molecules and the distance separating them [2,27,70]. 

In the Lennard-Jones model the sum of the attractive (London or dispersion forces) and repulsive forces and 
potential energy function for spherical molecules is given by 

F = K R r~ 13 + K L r ~ 7 (27) 

t/(r) = [hjf] 1 " -12 — [ — r_ 6 (28) 

In which r is the molecular distance. The Lenard-Jones potential function is commonly written in the form of 

«w = 4 £ {[yf 2 -[yfj (29) 

e is the maximum energy of attraction is one value of r for u(r) = 0. 

In their assumptions, the polymer molecules are as the gas molecules in the L-J model with a potential energy that 
varies with the distance between them. Solvent and particles also increases the intermolecular separation by 
occupying the interstitial space between the molecules. The effect on the free volume would be proportional to the 
solubility. Both of these effects act to increase the potential energy of the system. These effects has quantified into a 
coherent expression [2,27,70]. 

The molar volume of a polymer consists of four terms: a reference volume at a reference temperature and pressure, 
volume changes due to temperatures difference from the reference, any internal gas pressure and any additives in 
solution. 


V(T, P, C) = [V 0 + a(T - r 0 )] + ( 1 + + (Volume of additive ) (30) 

Where V 0 is the molar volume of the polymer repeat unit at a reference temperature, a is the coefficient of thermal 
expansion of the polymer, T 0 is the reference temperature, T is the temperature, P sat is the saturation pressure of the 
gas in the polymer , B(T) is the bulk modulus of the polymer at temperature T. 

Polymer molecules can be modeled as points with spherically symmetric potential fields. 


The volumetric molar mass is related to the distance between polymer molecules: 

a = 2 (— V w ) /s 
UttN w J 

R(T,P,C ) = l[- 

By replacing R in the potential energy relation we come to: 


3 >V3 


U(R ) = U 0 


C R %T-K Ro /r )1 


(31) 

(32) 


(33) 


Where Uo is the equilibrium potential energy. Thus, changes in the potential energy of the polymeric molecules 
originating from the presence of gas and additives can be calculated. Potential energy relative to changes in 
pressure and concentration of additives show a stronger and weaker dependency respectively and increase in the 
distance between polymer chains enhances the variations in the potential energy [26,48] . 


Colton and Suh included the potential energy in their calculations of nucleation energy but did not reach 
acceptable results which they assumed to originate from the incorporation of J-L potential energy model. The free 
volume model was proposed as the first methods for the modification of the thermoplastic foam version of the 
classical nucleation theory for resolving the mentioned inconsistency, although there have been criticizes on the 
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model and its assumptions [2, 27,70]. 

In another investigation, an attempt was made to reduce the errors of the classical nucleation theory for the 
prediction of energy and rate of nucleation in polymeric foams by Han and Han [40]. 

Which the effect of polymer chain and super saturation on nucleation energy can be expressed by the following 
relation: 

A G t = K B T[Lncp 1 + (p 2 + XV%\ (34) 

AG s = nK B TLn(j|) (35) 

Where and are the volume fractions of solvent and solute, respectively, x the interaction parameter, n is the 
number of solvent molecules in a critical bubble, PG is and pressure of the gas inside the nucluous, PV is the 
equilibrium vapour pressure and is the super saturation level. 

Thus, 


AG£ = AG* -AG t -AG S (36) 

Where, AG* is the critical energy of the classical nucleation theory. 

Finally they have calculated the density of cells theoretically and evaluate their consistency with experimental 
results. The theoretical results have higher values in comparance to experimental results. The reason for this 
inconsistency is recognized as the coalescence of nucleuses. Despite the unpleasantness of this outcome, usually the 
results gained from the classical nucleation theory exhibit smaller values of cell density and nucleation rate 
compared to the experimental results which stands as a proof for the effectiveness of the performed modifications 
[40,55]. 


Some of the studies for the evaluation of the parameters of the classical nucleation energy relation with more 
precise methods will be reviewed in this section. 

The variants of the nucleation energy relation are the difference between the interfacial tension of the polymer and 
the solvent nucleus and also the difference between the inside and outside pressure of a nucleus which are 
respectively equivalent to the difference between the macroscopic interfacial tension of a smooth surface between 
solvent and polymer and also the pressure difference of the saturation stage and atmosphere in foaming process. 
The macroscopic interfacial tension is strongly dependent on factors such as temperature, pressure and weight 
fraction of solvent in polymer. In order to estimate the solvent weight fraction (usually a gas), one has to calculate 
its density at a specific thermodynamic condition with the use of an appropriate equation of state [39,98-100]. 

Satish and his coworkers have developed polymethylmethacrylate foams in a batch process with the incorporation 
of carbon dioxide. They traced the effects of pressure and temperature changes in the ranges of 0 to 35 MPa and 40 
to 80 °C on the amount of gas absorption, interfacial tension between gas and polymer, nucleation energy and 
critical radius, theoretically using the modified mean field lattice- gas theory (MFLG). They also managed to 
calculate the nucleation rate experimentally and theoretically through the use of homogeneous classical nucleation 
theory. A summary of their results will be reviewed presently [39,101]. 


The following relation can be used for the calculation of the Interfacial tension between the polymer and gas phase. 


y 1 / 4 = 
1 mix 


o YL 5Zl 

unixAi-1 jj. 


1/4 


(37) 


Where, Xj, q. and y. are respectively the molar fraction, molar density and interfacial tension between the 
components of a mixture. Knowing that the interfacial tension of a supercritical fluid is zero, by replacing molar 
density with mass density (g/cm3) we have: 


Y n 


Y 


(- 


polymer \^ polyme] 


(l w gas) 


(38) 
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Where w gas is the weight fraction of the absorbed gas in polymer at saturation. With the calculation of weight 
fraction of the absorbed gas and mass density of the mixture using the MFLG model or any other suitable equation 
of state for the polymer-solvent solution at high pressure, the mixture interfacial tension will be reached and the 
values of critical radius, critical energy and nucleation rate at various temperature and pressure condition will be 
found. The number of nuclei can be estimated using the relation below: 


N, 


total 


= jtjdt = j Po j 

J 0 J J Psnt J 


dP 

dP/dt 


(39) 


Where Psat and Po are the saturation and atmospheric pressure respectively. 


In addition, the number of cells can be assessed experimentally using the SEM results with the relation below: 

^ 3/ 2 (40) 


N 


f nM 2 \ 
f ={~) 


Where n, M and A are the number of cells in the SEM image of the foam, magnification and the image area. 


According to the experimental results, pressure increase and pressure drop rate and also temperature reduction 
result in an improved nucleation. Changes in carbon dioxide absorption in polymer with pressure and temperature, 
changes in interfacial tension of polymethylmethacrylate with temperature variations of its glass temperature with 
pressure were discussed and the results have a decent consistency with the experimental data. Comparison 
between the number of nucleus from relations 50 and 51 indicate an acceptable consistency between the theoretical 
and experimental results especially at pressures above 20 MPa [39]. 

Thus, the results gained from this study confirm the possible reduction of the existing errors between the 
experimental and theoretical results in some experimental conditions with the precise calculation of the 
experimental parameters of the classical nucleation theory with a suitable equation of state. 

In a study on foaming of polystyrene and polylactide with carbon dioxide at various temperatures and pressures, 
the classical nucleation theory and the Non-Random Hydrogen-Bonding (NRHB) were used for the prediction of 
results [87]. 


The NRHB state equation is frequently used to investigate the thermodynamic behavior of gas and liquid phases, 
supersaturation state, polymeric solutions and also for the evaluation of gas absorption or supersaturated fluids in 
polymers and its effect on plasticization [98, 99]. 

Based on the experimental data, increase in pressure and pressure drop rate and decrease in temperature results in 
an improved nucleation. In addition, with the used of the NRHB model it has become possible to assess and 
monitor the changes in carbon dioxide absorption in polymers, variations of interfacial tension of polystyrene, 
changes in its glass transition temperature with pressure and temperature and the final results have demonstrated 
a decent consistency with the experimental data. 

The changes in nucleation energy, cell density and critical nucleus size with pressure and temperature and have 
also been studied using the values obtained from the thermodynamic model and the classical nucleation theory. 
Qualititative analysis of changes in nucleation energy in regard to the size of critical nucleus and a quantitative 
analysis of cell density in relation to experimental data appropves the capability of the classical nucleation theory 
for accurate prediction of results both quantitatively and quantitatively in case of using correct values of data [98]. 

Thus, in some cases the experimental data for microcellular foams have good consistency with the classical 
nucleation theory and its modifications, although the inconsistency increases with the reduction of cell size [50]. 


Density Functional Theory 

In the unstable region of the phase transition diagrams, fluctuations of temperature, density or pressure are in a 
limit where the energy barrier is conquered and phase transition occurs, while the classical nucleation theory does 
not predict a reduction of the energy barrier around or inside the unstable region. In the classical nucleation theory 
the thermodynamic properties of nucleus are calculated based on the thermodynamic properties of the bulk and a 
distinctive boundary is presumed between the nucleus and the bulk, while no exact boundary can be predicted 
physically. In this theory the interfacial tension is considered between the bulk and the macroscopic nucleus and is 
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calculated by assuming a smooth and curveless boundary which is the source of error. Therefore, in order to 
overcome such a conspiracy, the density functional theory can be incorporated for the modeling of nucleation [42- 
49, 102-103]. 


In statistical mechanics, the macroscopic thermodynamic methods and methods based on microscopic simulation, 
the density functional theory is a new method for the study of substances. In this technique the molecular details of 
substances are used for the development of a rational calculation time and an acceptable physical sense [102-103]. 

In this theory, nucleation is considered as a local density fluctuating of molecules and with the elimination of the 
dividing surface between the new and old phases; the free energy of the system is calculated based on the 
fluctuations of density or component compositions. 

If p(r) is regarded as the density of molecules at point (x,y, z) with a position vector of r(x,y, z), then space 
inhomohenety or fluctuations of the molecular density of the system is expected which results in the fluctuations of 
quantitative properties. The free local Helmholtz energy per molecule f is a function of r, density and its 
derivations and the free energy of system in the two phase condition is expressed as: 


f 2 (q) = J v f(r)p(r)dr ( 41 ) 

Which o(r)dr is the number of molecules in an infinitesimal differential volume of dr = dx dy dz around point r. 

For the free Gibbs energy we have: 

G 2 = F 2 + PV (42) 


Equation (52) is a function of density function. When a function is dependent on another function, it is known as a 
functional. 


In the density functional theories, the work of nucleation w is a function of density function. 

G 2 {q} = J v [ f ( r ) Q<T ) + Q]dr ( 43 ) 

For the work of nucleation we have: 


w{q} = f v { [f(r )- p old ] Q(r ) + Q ] } dr 


(44) 


Relation (55) is the work of formation of density fluctuation with a profile p(r) in a single component system 
which there is no difference between the new phase and the bulk [29, 30]. 

The critical work of nucleation w* , is a special form of work for the formation of density fluctuation: 

w* = w{q*} (45) 


Where, 

(dw {q})q=q* = 0 (46) 


p*(r ) is obtained from solving equation (57). 


The free energy function f(r) has to be known for solving equations (57) and (56). The most common way to define 
this function is through the potential molecular functions. 


Cahn and Hilliard used the density functional theory for the investigation of nucleation in a two component fluid 
for the first time in which the free energy is calculated based on the regular solution. Their results for critical 
nucleus were in contrast with the results obtained from the classical nucleation theory which by approaching the 
unstable region, nucleus size is increased while nucleation energy is reduced. 


This theory has been used in many cases for the evaluation of nucleation. In some studies, the effects of size and 
surface curvature of nucleus on interfacial tension has been calculated using this theory and the achieved results 
have been used in the classical nucleation theory. Additionally, the value of the nucleation energy calculated from 
this equation is found to be less than that of the classical nucleation theory. The theories based on this standpoint 
have also been used to assess the nucleation energy in polymeric foams. Since these theories are based on 
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molecular models, thus their results are only sufficient for investigating the effects of various parameters and 
insight into nucleation process and methods of promoting its efficiency. However, these results are not comparable 
to the experimental results. Because, first, the materials studied experimentally are relatively complex and may not 
be describable by simple spherical potentials. Second, nucleation rates are exponentially sensitive to almost 
everything in the calculation; a small error in calculating the surface tension from first principles could translate 
into a huge error in the predicted nucleation rate [3, 42-49, 104-105]. 


Self Consistent Field Theory (SCFT) 


The similar sizes of polymer chains with high molecular weight has manifested the role of global polymer behavior 
over the effect of atomic details for characterization and development of common theories for modeling of 
polymeric systems. Another feature of polymeric systems is the many number of configurations which results in 
contact between a polymer molecule and many other molecules and reduce the effect of intermolecular forces. 
These characteristics have made the incorporation of mean field theories possible for modeling and description of 
polymeric systems [103,106]. 

Self consistent field theory, the most applicable mean field theory with fundamental and reliable assumptions 
commonly is used for the study of polymeric blends, polymer absorption at surfaces, microphase separation in 
block copolymers and phase separation in blends and polymeric solutions. This theory is an expression of the 
changes in free energy of polymeric solutions with variations of composition which is regarded to be in 
competence with the density functional theory for the calculation of nucleation energy in polymeric foams 
[103,106-110], 


The free energies of a fluid bubble in a polymeric solution and single phase polymeric solution can be calculated 
using the SCFT with the following relations: 

NF 


T = 


pokgTV 


= ~ (1 - 0s)ln [n5y - J ln O + y dr [XN VOcMr) “ (47) 

w s (r)cf> s (r) - w p (r)cf> p (r)] 

= (l-0 s )ln(l-0 s )+^lnQ+ X N(l-0 s )0 s (48) 


T = 


e 0 k B Tv 


Where the cf> (r) and cf> (r) are the local volume fractions of polymer and solvent respectively, (1 — 0 S ) and 0 S are 
the global volume fraction of polymer and solvent respectively and Q p and Q s are respectively the partition 
functions of a single molecule of polymer and solvent. The volume ratio of a solvent molecule to a polymer 
molecule is characterized by a, the volume of a polymer segment is represented by qr 1 and the polymerization 
degree is accounted for by N. 


The Florry Huggins parameter % expresses the separation degree, chemistry of solvent molecules and segments of 
polymer and is inversely proportional to temperature. w s p (r) is a function of interaction between molecules of 
solvent and segments of polymer with other solvent molecules and other polymeric segments respectively, and 
also the type of mean field theory of SCF. 

The derivations of equation (57) in relation to the mentioned functions create a set of nonlinear equations and are 
solved using an appropriate equation of state with regard to a polymeric system through the self consistent 
method from which the energy of a two phase system is obtained [3, 38]. Furthermore, the incorporation of this 
theory in various studies for the analysis of nucleation phenomenon in polymeric solutions will be reviewed. 

The principal reasons for the inconsistency between the results of classical nucleation theory and experiments are 
neglecting the effect of cell size reduction and surface curvature on the interfacial tension and the error occurring in 
the calculation of volume energy difference of nucleus and bulk (pressure difference). 

The SCF theory is a microscopic model which includes the degrees of molecular freedom and is a suitable 
alternative for the classical nucleation theory. Therefore, in most studies the results of these two theories are 
compared and the insufficiency of the classical nucleation theory for a close prediction of results due to underlying 
microscopic reasons is discussed. Furthermore, the results of this theory are not usually comparable to the 
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experimental data and can merely be used for a qualitative analysis and evaluation of the trends of the effective 
parameters. Actually, SCFT provides an upper bound on the cell density, with experiment expected to find lower 
values of cell density. If one is interested in nanocellular foams for example, one can use SCFT to avoid spending 
too much time experimentally searching points of the parameter space for which SCFT maximal cell density 
predictions are low and concentrate on areas where the cell density has the possibility of being high[3, 38]. 

In one research, the interfacial tension and the volumetric energy difference are calculated using SCFT and 
subsequently the obtained figures are put into the energy relation of the classical nucleation theory. The 
application of the classical nucleation theory means accepting of the initial assumptions, i.e. a distinctive dividing 
surface between the nucleus and the bulk. The nucleation energy can also be calculated directly from solving 
equation (58) which represents the free energy of a solvent nucleus in a polymeric matrix. According to the figure 
(7), the classical nucleation theory overestimates the critical radius and the nucleation energy. 



FIG7. COMPARING OF CHANGE OF NUCLEATION ENERGY WITH NUCLEUS RADIUS IN CNT AND SCFT 

Actually the application of the SCF theory presents a possibility for introducing the effect of curvature of interface 
in scale of polymer molecules on the interfacial tension and nucleation energy. Increase of polymer molecules 
conformation (enhanced entropy) and reduction of internal energy reduce interfacial tension and nucleation 
energy [38]. 

In another research, the SCF theory has been used for the estimation of critical nucleation energy, critical size and 
density of cells. The principal purpose of this study was the achievement of an experimental method for the 
manufacture of a foam structure with high quality with unique characteristics such as smallest possible cell size, 
maximum interfacial area and nucleus with no polymeric chains. The effects of changes in the amount of the 
experimental parameters on the foam structure have been studied in this work. Owing to the fact that the SCF 
theory is based on the equilibrium statistical mechanics and does not generate any information on the kinetics of 
nucleation, it is therefore not a suitable method for the prediction of theoretical results and their comparison to 
experimental data. It can merely assess the conditions for the development of foam with an ideal structure from 
theoretical and qualitative point of views. In addition, the resultant information from the SCF theory are applicable 
for molecular analysis, deep understanding of nucleation phenomenon, identification of experimental methods and 
experimental conditions and reducing the number of experiments[3]. 

Conclusion 

Based on literature the physical foaming of thermoplastics consists of first forming a saturated mixture of solvent 
and the polymer and then by changing the free energy of the mixture by a pressure drop or temperature rise, 
producing a two phase mixture consisting of solvent molecules in gas state surrounded by a polymeric layer. In 
this process the change in free energy of the states, after and before the spherical gas nucleus formation in polymer 
bulk is calculated using nucleation theories. 
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The main theory that has been used for this process is the classical nucleation theory. The controlling terms in 
using this theory for thermoplastic foaming process are interfacial tension and pressure difference between gas 
nucleus and polymer layer. Inaccurate and/or incorrect calculation of these terms, like neglecting the effect of 
polymer chains on the free energy of the system or disregarding the effects of time and pressure drop rate, are the 
main reasons for the inconsistency between the results of this theory and the experimental data. CNT is commonly 
used for prediction of thermoplastic foaming nucleation rate and its energy terms are mainly calculated using a 
suitable equation of state. 

In the density function theory, the assumption of the existence of a dividing surface between the nucleus and the 
initial phase and the related errors are eliminated. In this theory, nucleation is considered as a local fluctuation 
density and there is no boundary between new and old phases. In addition, the free energy of the system is 
calculated based on the fluctuations of density of the components. This theory is used in many cases for the 
evaluation of nucleation phenomenon in phase transitions. In some investigations the effects of size and surface 
curvature of nucleus and experimental conditions on interfacial tension have been assessed using this theory and 
the obtained results have been incorporated in the classical nucleation theory. Additionally, the nucleation energy 
calculated directly using this theory appears to be lower than its equivalent calculated using the classical 
nucleation theory which can decrease the inconsistency between the results of this theory and the experimental 
data. It should be stated that, this theory has not yet been evaluated with reference to experimental results and it is 
merely used for the study of nucleation in molecular scales. This theory is been used recently in a preliminary work 
for prediction of nucleation in thermoplastic foam. 

The variations of free energy of a polymeric solution with changes of components composition is expressed 
through the mean field theory of SCF. This theory has shown the potential to incorporate the effects of interfacial 
surface curvature in the scale of polymer molecules and it has been used to assess the critical energy of nucleation, 
critical size and density of cells. The SCFT is not considered as an appropriate method for the prediction of 
theoretical results and its comparison to experimental data could merely be used for the evaluation of the 
conditions for obtaining foams with ideal structures from the theoretical and qualitative point of views. 

In computations using DFT and SCFT methods for homopolymer foams, the size of the formed nuclei in the 
nucleation stage is estimated to be less than 10 nm[3,38,104-105], while the results from the electron microscopic 
images of homopolymer foams indicate cell sizes of about 1 pm or more. If we accept the computational results of 
mentioned theories and assume that the nuclei formed in the nucleation stage in homopolymer foams are less than 
10 nm, the produced foams must have certainly passed the growth stage as well as the nucleation stage. Therefore, 
the available equipments do not have the capability to record the nucleation stage in homopolymer foams. Also the 
nucleation rate deduced from the classic nucleation theory are completely incorrect, because its content is even less 
than the reported experimental rate from the SEM images. 

It can be presumed that nucleation has taken place completely in homopolymer foams and the next set of events 
after the nucleation stage and during the growth stage determine the structure of the foam. In fact, solvent 
molecules are dispersed between polymer chains and have a much higher mobility. With the beginning of pressure 
reduction, in solvent molecules phase transition will take place from the super critical state with higher density to a 
gas state with lower density through nucleation mechanism and the initial nuclei will be formed with sizes less 
than 10 nm. 

Since first order phase transition of solvent with nucleation mechanism is the controlling factor for the nucleation 
stage, the surface tension between the gas and super critical state of the solvent molecules is negligible, resulting in 
a very small nucleation energy, thus, the nucleation efficiency and density can be considered very close or even the 
same in foams with equal solvent absorption. After nucleation, the shaped nuclei tend to grow and escape from the 
bulk to decrease chemical potential difference. 

The nuclei are placed in the free volume between polymer chains and in order to escape and grow they need to 
overcome the chain movement energy (chain elastic energy) and they make their way out through the movement 
of entanglements and chain slippage. Therefore there is a conflict between nuclei growth energy and elastic energy 
of chains where the domination of each of these energies determines the size and structure of the cells. Foam 
structure will be stabilized in case of temperature reduction in foam to less than the glass transition temperature of 
gas and polymer mixture. 
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Thus it seems that a model with the ability to predict the conflict between the two mentioned energies for the 
prediction of final cell density will be best suited. 
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